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CN ' Abstract 

We suggest modified bi-level approach for finding the best stacking sequence of laminated composite struc- 
tures subject to mechanical, blending and manufacturing constraints. We propose to use both the number 
f ) ■ of plies laid up at predefined angles and lamination parameters as independent variables at outer (global) 

stage of bi-level scheme aimed to satisfy buckling, strain and percentage constraints. Our formulation allows 
precise definition of the feasible region of lamination parameters and greatly facilitates the solution of inner 
level problem of finding the optimal stacking sequence. 



Introduction 

Due to the excellent mechanical properties laminated composite materials are widely used nowadays in various 



> 

pg industries (see, e.g., Ref. pQ for review). The common problem is to develop composite structures of minimal 

f^. , weight subject to mechanical, blending and manufacturing (technological) constraints, the degrees of freedom 
being the number of plies laid- up at a particular orientation (angle). Therefore, the industry faces with extremely 
difficult high dimensional mixed-integer optimization problem, in which both total plies percentage and the 
concrete stacking sequence of plies are to be found in order to get optimal design satisfying all the imposed 
constraints. 

It is important to note that substantial part of constraint functions are stacking sequence independent. 
Indeed, the weight of laminated composite structure is simply proportional to the total number of plies. On 
the other hand, imposed constraints could be naturally divided into two classes: i) "universal" mechanical 
loads written in terms of various stiffness tensors; ii) "non-universal" manufacturing constraints, which vary 
for different applications and dictate, for instance, particular laminate blending rules (see Ref. [5] for details 
and further examples). Let us consider the most important mechanical constraints, which according to classical 
laminate plate theory [2] (for review see, e.g., Ref [3] and references therein) are formulated in terms of in-plane 
A, coupling B and out-of-plane D stiffness tensors 

(1) 

linearly relating mid-plane strains e and plate curvatures k to vectors of running loads N and out-of-plane 
moments M . To simplify analysis we follow well-established industry requirements according to which only 
the case of symmetric laminates with finite (usually small) number of possible ply orientation angles 9{ G O = 
{0i, ..., OatqIi are to be considered. Then B tensor vanishes identically while in-plane and out-of-plane stiffness 
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matrices can be expressed as linear functions of four-dimensional lamination parameters £f, £P , i = 1, 2, 3, 4 

e 4 = /^ Ef =1 C? } A. 

?> = a^ELcT'^Az (2) 

Cj 0) = [cos(2ft), C08(4fii), sm{26 l ), sin(4^) ] T . 

Here N is the half of total number of plies each of thickness Az, with h = NAz being the panel half-thickness, Oi 
denotes orientation of i-th ply, z% stands for the distance between mid-planes of laminate panel and i-th ply and 
we used vector notations instead of explicit indices (this nomenclature is followed below). It is apparent that 
in-plane mechanical loads are independent upon the stacking sequence 9 and are the functions of lamination 
parameters £ only, while for out-of-plane constraints stacking sequence dependence is totally hidden in £ D 
vector. Note that £ lamination parameters are functions of only the number of plies n% laid up at particular 
angle 6i (ply numbers n). Hence dependence upon £ A could equally be written as n dependence, but reverse 
statement does not hold in general, the function £ (n) is invertible only in special case Nq = 4. As far as 
buckling load constraint is concerned it brings essentially nothing new to the above picture. Indeed, buckling 
load factor is linear in D matrix coefficients (see, e.g., Ref. [3]) and hence depends on stacking sequence only 
through the £ lamination parameters. 

Unfortunately, it is hardly possible to perform similar analysis for "non-universal" constraints mentioned 
above. These are extremely case-to-case dependent, include specific manufacturing requirements for each par- 
ticular product and thus are to be taken as most general stacking sequence dependent constraints. Therefore, 
the composite weight optimization problem can be summarized as 

min J2i n i 

P(n) > "universal" ; ply number dependent 

(3) 
s t ) L(£ D ) > "universal" ; lamination parameters dependent 

S(8) > "non-universal" ; stacking sequence dependent 

Suitable solution methodology crucially depends upon the problem size. It is true that for thin laminates 
consisting of only a few plies direct methods are applied well and outperform alternative approaches. However, 
in engineering applications ply number N for each plate might be of order one hundred, not mentioning the 
need to consider many different plates glued together. Then direct treatment becomes extremely inefficient and 
alternative methods need to be developed. Nowadays the most widely used approach is to apply approximation 
techniques (see Refs. [9] for review), in which "universal" constraints P, L are treated rigorously at outer level, 
while the last set of "non-universal" S'-constraints are accounted for only approximately at inner stage, iterating 
the whole procedure in case of large discrepancies between outer and inner levels. The possibility of problem 
decomposition into outer/inner levels arises from crucial observation that none of "universal" constraints depend 
on stacking sequence directly. Thus at outer iteration we can formally consider £ , £ D as independent variables 
and solve the auxiliary taslo 

[ P(ft) > 
min > Hi s.t. < (4) 

r 1 up) > o 

to get "optimal" lamination parameters £jf, £P , while at inner iteration we have to ensure that £, , £P are 
indeed realizable in term of particular stacking sequence 

min, |e - eVll 2 + II? ~ C°(0)\ 2 

(5) 
s.t. S(6) > 



1 Here we're slightly sloppy in notations and use ft instead of £ A in view of above mentioned correspondence between these 
quantities. 



Successful examples of the above methodology include bi-level composite optimization and lamination pa- 
rameters approach 8. , which differ only in formulation of outer level problem. In bi-level treatment one assumes 
high homogeneity of laminated composite and then derives simple proportionality of A and D matrices thus 
eliminating £ parameters from the list of design variables. Then the outer level problem becomes relatively 
easy to solve. However, these simplifications are not coming for free: once stacking sequence dependence is 
abandoned "optimal" ply numbers n* are not obliged to respect i-type (e.g., buckling) constraints, which is the 
prime disadvantage of bi-level optimization. Contrary to that, in lamination parameters approach one keeps the 
stacking sequence dependence explicit at outer level and considers lamination parameters £ , £° as independent 
design variables. It is crucial that no additional assumptions are introduced here, moreover, there is no need to 
check strain or buckling constraints as long as inner level optimization matches the optimal £* , £P values coming 
from outer iteration. However, the acute problem is to define the feasible region of lamination parameters (see 
Ref. 5 and references therein) so that the inner problem might be solved successfully. Unfortunately, there are 
only a few rigorous results on what is the feasible domain of £ A , £ D variables. Although it is known [10] that 
feasible domain of lamination parameters is convex, up to our knowledge no complete explicit equations tightly 
bounding allowed £ A , £ D values are available. 

Our approach to composite materials optimization lies essentially in-between the above bi-level and lam- 
ination parameters methods. Namely, we suggest to retain both £>-type lamination parameters £ and ply 
numbers n explicitly at outer optimization level, so that the formulation (j4j remains valid verbosely. As far 
as feasible region of £-° values is concerned, we show that at any fixed ply numbers feasible domain O^(n) of 
realizable lamination parameters £° is convex polyhedral body with Nq\ vertices, each of which corresponds 
to a particular "extreme" stacking sequences compatible of given ply numbers. Therefore, our approach is 
essentially the equation Q supplemented with feasibility requirement 

i° 6 ««(**) (6) 

and explicit description of the feasible region f^(n) as the convex hull of Nq\ points or equivalently as a set of 
linear constraints: 

n c (n) = {£:M D < b}. (7) 

Feasible Region of Lamination Parameters 

Derivation of the explicit form of lamination parameters feasible region n^(ft) at fixed ply numbers n = 

{ni, ..., 7ijv e } is similar in spirit to what had been done in seminal paper 7. Namely, for given vector A 

consider the problem 

0, = argmax (g° • A) , 9 = {9 1; ..., G Ns ] 

^e (8) 

which is an elementary step of constructing outer polyhedral approximation to f2g(n), [11]. Here 0^ denotes 
allowed ply orientation angles, S(x) stands for Kronecker delta-function and component-wise representation of 
£ is given in ([2]). We assert that stacking sequences, which deliver the maximum to ([8]) for various inputs A, 
are such that plies of the same orientation are adjacent 

0* w = {e fcl ,...,e fcl ,e fc2 ,...,e fc2 ,---e fcjVe ,...,e fcjVe }, » = i,...,N e \ (9) 

V v ' V v ' V v , 

™ fc l ™ fe 2 ™fcjv e 

Indeed, let us consider z-th term in the objective function (|5J|: 

z 1 ■ [A cos(26i) + Ai cos(40i) + A 2 sm(29i) + A 3 sin(4^)] (10) 

Suppose that it is maximized at a particular "best" ply orientation angle O^.. If there would be no constraints 
in the problem, solution must have all plies aligned at the same angle as well (these solutions constitute the 
vertices of famous Miki tetrahedron, Ref. [7])- In our case imposed constraints limit the maximal allowed 
number of plies with "best" orientation and actual task is to place these plies optimally in the stack. Due to the 



positivity of measure factor z 1 , which is maximal at composite skin layer and vanishes at laminate mid-plane, 
it is evident that the preferable position for ply with 0^ orientation is near the composite skin. Repeating the 
same arguments for all n^ plies laid-up at Ok we conclude that they occupy continuous stack located near the 
boundary layer. For the remaining ply orientations the same reasoning applies verbosely, the only difference is 
that boundary layer is now located at smaller z. 

It follows immediately from the above that all realizable values of lamination parameters £ D are located 
within the convex hull of "extreme" stacking sequences ©: fl^(n) 6 Conv(#* ). Confronting this with known 
concavity of lamination parameters feasible domain |10j we conclude that £ feasible region at fixed ply numbers 
must coincide with Conv(#* ) 

Q s (n) = Cohv(0* (m) ) . (11) 

This is the prime theoretical result of our paper, which opens the possibility to construct efficient bi-level 
optimization scheme for laminated composite optimization problem ([3]). Indeed, Eq. (JTTJ) might be used to 
ensure that outer level iterations always produce realizable optimal lamination parameters so that no difficulties 
can arise at inner level. Moreover, the number of vertices in corresponding convex polytope is relatively small 
and admits its efficient description in terms of linear inequalities. Indeed, in the most practically important case 
we have iVe = 4 and hence only Nq\ = 24 vertices, for which convex hull construction rises no computational 
issues whatsoever. 
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